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Abstract 

o 

Cn| a periodic Ising model has interactions which are invariant under translations of a 

full-rank sublattice £, oi 7?. We prove an exact, quantitative characterization of the 
Oh critical temperature, defined as the supremum of temperatures for which the spon- 

taneous magnetization is strictly positive. For the ferromagnetic model, the critical 
temperature is the solution of a certain algebraic equation, resulting from the condi- 
tion that the spectral curve of the corresponding dimer model on the Fisher graph has 
a real zero on the unit torus. With our technique we provide a simple proof for the 

'" ^ ^ exponential decay of spin-spin correlations above the critical temperature, as well as 

I the exponential decay of the edge-edge correlations for all non-critical edge weights of 

the corresponding dimer model on periodic Fisher graphs. 

B 

L_j 1 Introduction 

Phase transition is a central topic of statistical mechanics. A system undergoes a phase 
transition whenever, for some value of the temperature and other relevant thermodynamic 
parameters, two or more phases can coexist in equilibrium. Different properties of the 
pure phases manifest themselves also as discontinuities in certain observables as a function 
of the appropriate thermodynamic variables, e.g. discontinuity of the magnetization as a 
function of the magnetic field in a ferromagnet. Of special interest is the temperature at 
"O which the phase transition occurs, that is, the critical temperature. 

Lebowitz and Martin-Lof [22] defined the critical temperature, T^, for Ising ferromagnets 
^ ^ as the supremum of temperatures such that the spontaneous magnetization is strictly pos- 

^ itive. Lebowitz [20j identified Tc with the self-dual point for an isotropic, two-dimensional 

^ square grid Ising model. Based on various correlation inequalities and differential in- 

equalities, Aizenman, Bursky and Fernandez [2] characterize the phase transition by the 
exponential decay of two-point spin correlations above Tc, for ferromagnetic Ising models 
with dimension D >2. 

For the two-dimensional Ising model, another approach is to apply the Fisher correspon- 
dence [7j, which is a measure-preserving bijection between the even spin-spin correlation 
functions of the Ising model on a graph G, and the edge probabilities of a dimer model on 
a decorated graph, the Fisher graph. Since then, dimer techniques have been a powerful 
tool in solving the two-dimensional Ising model, see for example the paper of Kasteleyn 
[H], and the book of McCoy and Wu [25]. However, due to the complexity related to 
large matrices, the two-dimensional Ising model has acquired a notorious reputation for 
difficulty. On the mathematical side, Kenyon |16l 117^ [T8] has proven numerous spectacular 
results about the dimer model on bipartite graphs in recent years. Those results make it 
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possible to look at the two-dimensional Ising model from a new perspective. In our paper, 
we follow the dimer approach to the Ising model, and explicitly characterize the critical 
temperature, defined as in [22] and [2| , with the zero of an algebraic equation. 

Theorem 1.1. Let Tc be the critical temperature of the ferromagnetic, two-dimensional 
periodic Ising model, defined by the largest temperature such that the spontaneous magne- 
tization is strictly positive, then < Tc < oo. Tc is determined by the condition that the 
spectral curve of the corresponding dimer model on the Fisher graph has a real zero on . 

Theorem 1.1 allows one to compute the critical temperature of arbitrary periodic Ising 
ferromagnets accurately, by solving an algebraic equation, see Example 6.14. The proof of 
Theorem 1.1 consists of 3 steps: 

Step 1: Applying Lebowitz's technique [21] and the FKG inequality, we prove that 
the weak limit of the even spin-spin correlation functions is independent of the boundary 
conditions as the size of the graph goes to infinity. For the uniqueness theorem about the 
Gibbs measure of dimer models on a more general class of non-bipartite graphs, we refer 
to Corollary B.l. 

Step 2: Using an n x n torus to approximate the infinite bi-periodic graph, we express 
the two-point spin correlation functions as the determinant of certain block Toeplitz matri- 
ces. For basics of Toeplitz and Hankel matrices, we refer to the appendix. We prove that 
the determinant of the symbol of such block Toeplitz matrices is identically 1. Using the 
FKG inequality, we prove that there is a unique Tc,p, defined to be the lowest temperature 
such that the limit of the two-point spin correlation function is 0. The analytic property 
for the limit of two-point spin correlation function in case the entries of the block Toeplitz 
matrices are analytic with respect to the reciprocal temperature follows from Widom's the- 
orem [28] and a result (Lemma 4.6) from operator analysis [lOj. Hence Tc^p has to satisfy 
the condition that the spectral curve has a real zero on T^. 

Step 3: Characterization of Tc and the condition that the spectral curve has a real 
zero on (critical dimer weights). Since the two-point spin correlation functions are 
independent of different translation invariant Gibbs measures, we derive Tc^p = Tc, and at 
the critical temperature of the Ising model, the corresponding dimer model is also critical. 
We also prove that as f3 (reciprocal temperature) increases from to oo, there is a unique 
/3o (0 < /3o < oo), such that the spectral curve has a real zero on T^. The uniqueness of 
the critical dimer weights implies its identification with the criticality of the Ising model. 

Here is a description of the structure of the paper. Section 2 proves the uniqueness of 
even spin-spin correlation functions for periodic ferromagnetic 2D Ising models. In section 
3 we discuss the correspondence between Ising models on a square grid and dimer models 
on the Fisher graph. We also prove that the spectral curve is invariant under the duality 
transformation. In Section 4, we complete Step 2 and Step 3 to prove the main theorem. 
The block Toeplitz and the transfer matrix technique also provides a simple proof of ex- 
ponential decay of the spin-spin correlations at high temperature for Ising models with 
interactions satisfying certain symmetry properties, see Section 5. For a proof of exponen- 
tial decay for general interactions, we refer to ^ [1] . In the appendix, we discuss several 
basic facts of Toeplitz and Hankel matrices for the readers' reference. 
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2 Translation Invariant Gibbs Measure 



2.1 F.K.G Inequality for Periodic Interactions 

We consider a Ising spin system with a ferromagnetic pair interaction in a finite box A of 
Z^, embedded on the whole plane, i.e, at each point p of the lattice there is a spin dp = ±1, 
and the conditional probability of a spin configuration in the box A given a configuration 
outside it is proportional to 

where Jpq > is the pair interaction. Assume Jpq has period {m,n), i.e, for any i,j £ Z, 
Jpq = Jp+{im,jn),q+{im,jn)- ^ is the Uniform external magnetic field, and /3 is the reciprocal 
temperature. A boundary condition for the box A, is specified by giving a probability 
distribution 6a for the configurations outside A. Let T be the set of all Ising spin con- 
figurations on A. Under the partial order relation the partial order set F is a 
lattice, since any two configurations x and y in F have a least upper bound x V y and a 
greatest lower bound x A y. Moreover, F is a finite distributive lattice, see [8] for a 
definition. We have the following theorem: 

Theorem 2.1. Let f be an increasing function with respect to the partial order relation 
on T, fii,H2 be two probability measures, if all the following conditions are satisfied 

• the boundary condition 6i > S2; 

• the external magnetic field hi > /i2; 

• the pair interaction J^^ > J^^; 
Then the expected values of f satisfy 

Proof. Since F is a finite distributive lattice, according to |12j . if any spin configuration 
A,B satisfy 

/Xi(^VB)/i2(^A5) >^i(A)^2(5), (2) 

then 



for any increasing functions / on F. 

Let us prove Condition O) for 5i > 82, while hi = /12 and J^^ = J^^. If A does not 
have boundary configuration oi or B does not have boundary configuration 62, Condition 
Q is satisfied automatically since the right hand side is zero while the left hand side 
is nonnegative. Hence to prove Condition ([2]), it suffices to prove that for any A with 
boundary configuration 5i and B with boundary configuration 82, 

fi{AV B)fi{AAB) > fj.{A)fj,{B), (3) 

where fj is the probability measure on F with free boundary condition. Condition ^ is 
true if energy for all bonds Je > 0. Notice that Je may be different from edge to edge. It 
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is trivial to check Condition (|3]) for each single bond with 16 different combinations of spin 
configurations for A,B. If the inequality is true for the contribution of each edge in the 
probability, the inequality is true for the probability, since the probability is the product 
of contributions of each edge, and the partial order relation is defined bond by bond. 

The same method can be applied to check Condition Q spin by spin or bond by bond 
for the cases hi > /12 and J^^ > J^^. □ 

Remark. Condition (2) is called the F.K.G inequality. 



2.2 Free Energy 

As before, let A be a finite box of an infinite periodic square grid. The equilibrium state of 
the system in A is the probability distribution for configurations in A defined by ([T]) together 
with boundary condition 6a, or equivalently, by the family of correlation functions 

{(^A)h,A,bA = (Jl ^p)h,A,bA^ fo"^ ^ C A. 

An equilibrium state of the infinite system is defined to be a family of correlation functions 
obtained as the limit of correlation functions for a sequence of finite boxes with some 
boundary conditions 

{(^A)h,b = lim {aA)h,A,bA- 

A— >oo 

We say a sequence of boxes {An}^^i tends to infinity if 

1. each An is connected and A„ C A„_|_i; 

2. U„A„ = T; 

3. for any / > 0, let A' be the set of points in A with distance at most I from its 
boundary, and |A| denote the number of points in A, then 

|A^ I 

lim " = (Van Hove convergence) 

I I 



n— >oo 



The letter b indicates some arbitrary boundary condition, and b\ = " + " or 6a = " — " 
will indicate the boundary conditions defined by putting all spins outside A " + " or " — " . 
The statement {crA)h,b does not depend on the boundary condition is equivalent to the 
statement that the equilibrium is unique. 

For a finite system, the free energy is defined as follows 

F(/i,A,6A) = ^log(5]e-^AW), 

where E\{a) is defined in ([T]). We have the following theorem about the free energy 

Theorem 2.2. As A — )• 00, the limit of F{h,A, b\) exists and is independent of the bound- 
ary condition. Assume 

lim F(/i,A,6a) = F{h). 

A-¥oo 

Then F{h) is convex, and analytic for h ^ 
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Proof. First of all, for each fixed A and 6a, F{h, A, 6a) is a convex function of h. Assume 

= e-^^^"'^ > 0, 



peA 



Then 



> 0. 



The inequality follows from applying the Cauchy-Schwartz to {\/At}o- and {\/^A^B„}„. 
Equality holds only if all the B^r are equal, which is impossible, hence ^ ^^gf^'^^^ > 0, 
F{h, A, 6a) is convex. 

We claim that if limA-5>oo F{h, A, 6a) exists, it is independent of the boundary condition 
6a- Assume G{h, A) is the free energy for the free boundary conditions, namely 

G(/i,A) = — log^e^^^^p^'j6A,ip-q|=i-^p<j'^p'^9+EpeA'»'^p)^ 



then 



i^/? max Jpq < F{h, A, 6a) < G{h, A) + max 

|A| p,q \A\ p,q 



Under the assumption that A — )• oo in the sense of Van Hove, 



\dA\ 
lim ^—^ = 0. 

A^oo |A| 



Since Jpq has finite period, we have 



lim F(/i,A,6a)= lim G{h,A), 

A— !>oo A— s>oo 



which is independent of boundary conditions. 
It is convenient to introduce a new variable 



then 
where 

Q(z) = J;Q„z". (n = 0,l,...,|A|) (4) 

The coefficients Qn are the contribution to the partition function of the Ising lattice in zero 
external field from configurations with the number of " — " spins equal to n. The following 
lemma is proved by Lee and Yang |.23j 
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Lemma 2.3. Let x^p = 2:^^(0 7^ /3,a,/3 = l,2,...,n) he real numbers whose absolute 
values are less than or equal to 1. Let V = {1,2,-- - Divide the integers l,2,...,n 

into 2 groups a and h so that there are 7 integers in group a and [n — 7) integers in group 
b. Consider the product of all Xap where a belongs to group a and j3 belongs to group b. 
We shall denote by Qy the sum of all such products over all the ^|^^_|_^^, possible ways of 
dividing the n integers, in other words 



= Yi n n 

acV,\a\=-y a&a i3£V\a 



Consider the polynomial 



Q{z) = l + Qiz + ... + Q„z". 

Then all the roots of the equation Q = are on the unit circle. 

Apply the lemma to Q, we obtain that the roots of Q{z) = are on the unit circle. The 
distribution of roots of Q{z) = as A — 00 may be described by a measure dfi{6), so that 
\A\dfi{d) is the number of roots between e*^ and e*^^"'"'^^^ Since Q{z) has real coefficients, 
we have 

dfi{9) = dn{-e). 

Taking the logarithm of Q, and taking the limit as A — )• (X), we have 

f-2lT 

F{h) = /i/3+ / \og{z - e''^)d^l{e) 
Jo 

/"TT 

= hf3+ log{z^ -2zcos9 + l)dti{e). 
Jo 

The singularities of F{h) corresponds to zeros of z^ -2z + 1. Since z = e''^^^ > 0, the 
only possible singularity of F{h) happens at h = 0, z = 1. The convexity of F{h) follows 
from the convexity of F{h, A, 6a)- D 

2.3 High Temperature 

Define the average magnetization as follows 

. , . , 1 ^r- . l dFih,A,bA) 
m{h,KM) = (^L^P)'^A.A = dh • 

We have the following lemma about the average magnetization 

Lemma 2.4. As K ^ 00, the limit of m{h, A, b\) exists and is independent of boundary 
conditions. Moreover, 

lim m(/i,A,6A) = V/i 7^ 0. 

A^oo dh 

Proof. 1^ J2peA^p increasing function on Fa, and the F.K.G inequality says that 

m(/i,A, -) < A,6a) < A,+). (5) 
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For any A' 5 A, m(/i, A, +) can be obtained from m{h, A', +) by adding an infinite positive 
external magnetic field on A' \ A. By the F.K.G inequality we have 

m(/i,A',+) < m(/i, A,+). 

Hence m{h, A, +) is decreasing as A increases. Therefore limA_!.oo m^h, A, +) exists, denoted 
by m{h,+), for all h. Moreover, since \m{h, A,b\)\ < 1, according to the Dominated 
Convergence Theorem, we have 

f-h nh 

/ m{h,+)dh= lim / m{h, A, +)dh = F{h) - F{ho), 
J ho J ho 

where /i > /iq > or > /i > /iq- Since F{h) is analytic in h when h ^ 0, we have 

,^ , dF{h) 

Similar process shows that 

And the lemma follows from[5l □ 

Since F{h) is analytic in h when /i 7^ 0, we define 

* dF{h) 

m = hm . 

/i-s>o+ dh 

Applying F.K.G inequality and follow exactly the same process by Lebowitz and Martin-Lof 
[22], we have the following lemma: 

Lemma 2.5. When the external magnetic field is zero there is a unique equilibrium state 
for the periodic, ferromagnetic infinite system if and only if m* = 0. 

\X\ SpeA "^p increasing function on F, according to the F.K.G inequality, m{h, A, 6a) 
X^peA ^p)A,feA increasing function in Jij and h, hence if we fix Jij and h, it is an 

increasing function in /3. When /? = 0, we have a uniform distribution for all configurations, 

therefore 

m(/i. A, 6a)|/3=o = 0. 

As a result m{h, A, b\) > for any /3. As its limit m* is nonnegative and increasing in f3. 
The critical temperature f3c is uniquely defined by the conditions m*{f3) = for f3 < f3c 
and m*(T) > for /3 > /3c. Then we have the following theorem 

Theorem 2.6. When f3 < /3c, there exists a unique probability measure of ferromagnetic, 
periodic Ising spin system without external magnetic field. 
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2.4 Uniqueness of Spin-Spin Even Correlations 

This section is devoted to the proof of Theorem 2.7. The proof is divided into proving 
several lemmas, as specified below. 

Theorem 2.7. Given a bi-periodic, ferromagnetic Ising model with fixed reciprocal tem- 
perature (3, periodic interactions J, and external magnetic field h = 0, the spin-spin even 
correlation functions are unique under any translation invariant Gibbs measure. That is, 
fix a finite subset A of the vertices of the graph such that \A\ is even, for any translation 
invariant Gibbs measure /x, {cTAjfi is independent of n. 

Define 



PA 



n 



, 2 



Lemma 2.8. {pa)± = limA-s>oo(Pyl)A,± exists and is translation invariant. That is {pA+g)± - 
{pa)±, where g is a translation vector. 

Proof. The F.K.G inequality applies to / = because this is an increasing function 
in each spin variable. This implies that {pa)a',+ ^ {pa)a,+ , if A' 5 A, by the F.K.G 
inequality. Hence 

{pa)+ = lim {pa)a,+ 

A— s>oo 

exists and {pa)+ ^ {pa)a,+- The translation invariance of the {pa)± follows from the 
uniqueness of taking the limit A — oo. □ 

Lemma 2.9. 

{o'A)A,bA < {(^a)a,+ - 

Proof. Let k denote any edge connecting one vertex in A and one vertex outside A. Consider 



{<yA) 



z 

{'^a)' = J, . 



where 



2 = E 



a 



then 



(c^a) - (c^a) = (6) 

Y^{l — tA) ^ (JAe^^~J.'J6A0 Ji3{'^+Utj)a,aj+Y^^^Qi^{.Jk+tkJ',,yk 

" ~ 'ZZ' ' 

where ti = CTicr^. Here J'^ = iJk reflects the boundary conditions. It follows from the 
ferromagnetic condition that when 

Jk>\J'k\, ykedA, 
the right side of ([T]) is nonnegative. □ 
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Let A = BAC = BuC\BnChe the symmetric difference between the sets B and 
C, then tA = tstc = ±1, and 

l-tBtc = ±{tB-tc). (8) 

Substituting ([s]) on the right of ([t]) and going back to the a' variables, we obtain our basic 
inequahty for a a = ctbctc- 

{cJB<yc) - i^BCJcY > \{^B){acy - {^By{cJc)\- (9) 

Then absolute value sign on the right side of Q follows from the fact that the right side 
of ^ is always nonnegative given ferromagnetic interactions. Namely, Z, Z' are partition 
functions of Ising models, which are positive. For any give ti at all vertices, we consider 
an Ising model with interaction constants 

Jij = + titj) > 0, Vi ~ j. 

Jfc = + tkJ'k > 0, if Jk > I J'kl yk G OA. 

Since all interaction constants are nonnegative, we have a ferromagnetic Ising model, as a 
result, 

{aA)Z = (j^e^'~^'''j6A0 ■^'j(^+***j)'^''^j"+^fcs8A(-^fe+*fe4)'^ft > 

where Z is the partition function for the new Ising model. Moreover 1 — is always 
nonnegative, that is the reason that the right side of ([T]) is nonnegative, and ^ follows. 



Lemma 2.10. 



^l|m_^(PA)/3',+ = {pa) 13,+ , (10) 
^1™_(m)/3',- = {pa)p,-- (11) 



Proof. Since {pa)i3,+ < (pa) 



^lim^(pA)/3',+ < ^lim^(p^)/3',A,+ = {p. 



Letting A — t- oo, we have 
But F.K.G inequality implies that 



^l|m_^(PA>/3',+ < {pa)i3,- 



{pa)p,+ < {pa)i3',+ , 

for P < (3', so 

{pa)i3,+ < ^lmi^(pA)/3',+ . 



Then (10) is proved, and (11) can be proved analogously. □ 



Lemma 2.11. 



^ is continuous for all the ferromagnetic interactions. 
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Proof. First of all, notice that {aiaj)j^ = ((Tj(Tj)_ by symmetry. Let -§^y§^ 
right(left) derivative of F with respect to /3. The assertion will follow if we prove that 



be the 



dF 



lim 



(9^± A^oo |A| 



(12) 



where Ai is the quotient graph of the infinite periodic graph with respect to translation. 
Now we prove (12). Consider first the boundary condition +. For any e > 0, let be 
a box containing all the edges of Ai, such that {aiaj)\^^+ < {aiaj)+ + e, for all i ~ j, 
i,j £ Ai. Then for any translation vector g, and A such that A^ + 5 C A, we have for all 
i ~ e Ao, 



{aiaj)+ = {(Ti+g(7j+g)+ < (cJi+gCrj+g)A,+ 

Using the decomposition 



(13) 



where the second sum is over g, lim\. 



^ - n and 



Q imply 



that 



^ {aiaj)+ < liminf 

< lim sup 1^ ^(o-iO-j)A,+ < 7-^ Y + ^' 



|A| 



i,jeAi,i~j 



for any e, which means that 

1™ m E 



r dF{A,+) 1 
— , . , / \"t"'ilA + = hm — = 7- — r > [o'tCTj) 



and similarly for the boundary condition "— ". Since F is convex with respect to /3, the 

a/3 



discontinuous points of ^ is countably many /3's, at most. For any finite box A of a locally 



finite periodic graph 



«~i,j,ieA 



< maxdeg(f). 

I'SAi 



According to the Dominated Convergence Theorem, 

lim ^ ' df3 = lim 

A-5>oo ap A-5>oo 



1 dF{A,- 



dp 



-d(3 



= hm [F(A,/3i,+)-F(A,/3o,+)] = F(/3i)-F(/3o). 

A— s>oo 

The last equality follows from the fact that F is independent of boundary conditions. Let 
f3' ^ f] from the right side (/3' > /3). We have 



dF 



lim 



5/3+ a/3 



lim , . 

Ai 



i~i.«.ieAi 



So (12) is proved for the boundary condition +, and — is treated analogously. 



□ 
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Lemma 2.12. {aiGj)^ = {aiaj)j^ for any translation invariant limit probability measure /i 
and any i ^ j, if ^ is continuous. 

Proof. Let K be an arbitrary translation invariant interaction, define 



«~i.i.ieAi 



Then for any translation invariant probability measure fi 

F{J + K)>F{3) + {AK)f,. (14) 

That is (^k)/^ is a tangent vector of -F at J along the direction K. In other words, the 
expected value of Aj^ under any translation invariant probability measure corresponds to 

a tangent vector. For the proof of (14), see [19|. If ^'^^a^*'^^ is continuous, the tangent 

I I _ i=0 

vector is unique and its slope is equal to the corresponding derivative, which means the 

expected values of under any translation invariant probability measure are the same. 

Hence if ^ is continuous, we have 



E 



Given Lemma 2.9, under the assumption that Jij > 0, we have {aiaj)+ = {aiUj)^ for all 
translation invariant limit measure //. □ 

Lemma 2.13. Assume Jij > 0, if {(Jiaj) = {(Jiaj)' / 0, for all i ~ j, then {ue) = We)' , 
for all sets E containing an even number of sites. 

Proof. Any translation invariant Gibbs measure can be considered as the weak limit of 
Boltzman measures on finite graphs with given boundary conditions, as the size of the 
graph goes to infinity. Any boundary conditions can be transformed to " + " boundary 
conditions by changing the coupling constant on the edges incident to an boundary site. 
If the configuration at the boundary site is " + ", then the new coupling constant on 
incident edges is unchanged; otherwise the new coupling constant is —Je- Given Jij > 0, 
we have Jij > \Jij\, we derive that the inequality ^ is true. Since {asc^c) — {(^B(^C')' ^ 
\{(^B){cTcy - M'{ac)\, let A = BAC, we have 

(<Ta) - {ctaY > \{(^B){(^A(yB)' - (0-b)'(cJACTb)|. 

Then (cja) = {c^a)' , (o"_b) = (cfi)' / implies {ctactb) = {cactb)' ■ The result follows by 
induction. □ 

Proof of Theorem 2.7 From Lemma 2.13, it suffices to prove that (fTjCTj) 7^ and is 
independent of boundary conditions for all i ~ j. From Lemma 2.12, this is true if and 
only if ^ is continuous, and ^ is always continuous, given lemma 2.11.0 

The uniqueness of the spin-spin even correlation functions can lead to the uniqueness 
of translation invariant Gibbs measures of dimer models on a large class of non-bipartite 
graphs. Interested readers may look at the appendix for more about this topic. 
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3 Fisher Correspondence 



3.1 Dimer Model 

The Fisher graph we consider in this paper is a graph obtained from a honeycomb lattice 
by replacing each vertex by a triangle, as illustrated in Figure 2. There are three types 
of non-triangle edges on such a graph with different direction, namely o-type, 6-type, and 
c-type. We give weight 1 on all the a edges(corresponding to horizontal edges in Figure 1). 
The weights of 6-edges and c-edges are strictly less than 1 and assigned periodically with 
arbitrary period mxn. As we will see, the dimer model of such a Fisher graph corresponds 
to the ferromagnetic Ising model on a periodic square grid. 




Figure 1: Fisher Graph 

A perfect matching, or a dimer cover, of a graph is a collection of edges with the 
property that each vertex is incident to exactly one edge. For general results of dimer 
models, see the Appendix. 

A perfect matching, or a dimer cover, of a graph is a collection of edges with the 
property that each vertex is incident to exactly one edge. 

To a weighted finite graph G = (V, E, W), the weight W : E ^ M"*" is a function from 
the set of edges to positive real numbers. We define a probability measure, called the 
Boltzmann measure n with sample space the set of dimer covers. Namely, for a dimer 
cover D, 

f^{D) = ^ n ^(^)' 

where the product is over all edges present in D, and Z is a normalizing constant called 
the partition function, defined to be 

^ = E n ^(«)' 

D eeD 
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the sum over all dimer configurations of G. 

For planar Z^-periodic graphs, endowed with periodic weights on edges, the notion 
of Boltzman measure is replaced by that of Gibbs measure with the property that if 
the dimer configuration is fixed in an annular region, then the random dimer configuration 
inside and outside the annulus are independent, and the induced probability measure inside 
the annulus is the Boltzmann measure defined above. It is known that different Gibbs 
measures may be obtained as the infinite- volume, weak limits of Boltzmann measures with 
various boundary conditions. 

If we change the weight function W by multiplying the edge weights of all edges incident 
to a single vertex v by the same constant, the probability measure defined above does not 
change. So we define two weight functions W, W' to be gauge equivalent if one can be 
obtained from the other by a sequence of such multiplications. 

The key objects used to obtain explicit expressions for the dimer model are Kasteleyn 
matrices. They are weighted, oriented adjacency matrices of the graph G defined as 
follows. A clockwise odd orientation of G is an orientation of the edges such that for 
each face (except the infinite face) an odd number of edges pointing along it when traversed 
clockwise. For a planar graph, such an orientation always exists |14j . The Kasteleyn matrix 
corresponding to such a graph is a x |y(G)| skew-symmetric matrix K defined by 

{W{uv) if n ~ V and n — )■ v 
—W{uv) if u ~ V and u <^ v 
else. 

It is known [131 [T^ [271 E] that for a planar graph with a clock- wise odd orientation, the 
partition function of dimers satisfies 

Z = Vdet K. 

Now let G be a Z^-periodic planar graph. By this we mean G is embedded in the plane 
so that translations in 7? act by weight-preserving isomorphisms of G, i.e. isomorphisms 
which map each edge to an edge with the same weight. Let G„ be the quotient graph 
G/(nZ X nL). It is a finite graph on a torus. Let 7a;,n(7y,n) be a path in the dual graph 
of Gn winding once around the torus horizontally (vertically). Let EniEy) be the set of 
edges crossed by ^x,n{ly,n)- We give a crossing orientation for the toroidal graph Gn as 
follows. We orient all the edges of G„ clockwise odd except for those in Eh U Ey- This 
is possible since no other edges are crossing. Then we orient the edges of Eh clock-wise 
odd as if Ey did not exist. Again this is possible since G — Ey is planar. To complete the 
orientation, we also orient the edges of Ey clockwise odd as if Eh did not exist. 

For 9,T ^ {0, 1}, let Kn^ be the Kasteleyn matrix in which the weights of edges in 
Eh are multiplied by (—1)^, and those in Ey are multiplied by (—1)^. It is proved in [27] 
that the partition function Z„ of the graph G„ is 

Zn = ^|Pf«) + Pf«°) + Pf(^^') - Pf(i^^')|- 

Let Em = {ei = uivi, em = UmVm} be a subset of edges of G„. Kenyon [T6] proved 
that the probability of these edges occurring in a dimer configuration of Gn with respect 
to the Boltzmann measure P„ is 

Pnieu...,em) = nr=l ^(^.^0 ipf(^OO)^^ ^ Pf(K^°),,^ + Pf«^)i.^ - mK'n')E^J, 



13 



where E!f^ = V{Gn)\{ui,vi, ...,Um,Vm}, and {K^'^)e^ is the submatrix of K^"^ whose hnes 
and columns are indexed by E^. 

The asymptotic behavior of Zn when n is large is an interesting subject. One important 
concept is the partition function per fundamental domain, which is defined to be 

1 

lim (Zn)^ ■ 

n— >oo 

The logarithm of the partition function per fundamental domain is the free energy, 
namely 

F := lim ^logZ„. 

n— >oo n 

Let Ki be a Kasteleyn matrix for the graph Gi. Given any parameters z, w, we construct 
a matrix K{z,w) as follows. Let "fx,!, 7j/,i be the paths introduced as above. Multiply 
^i.{u,v)i the {u,v) entry of the matrix Ki, by z, if the crossing orientation on that edge 
is from u to v, otherwise multiply j-^ ^) by ^, and similarly for w on jy. Define the 
characteristic polynomial P{z,w) = det K{z,w). The spectral curve is defined to be 
the locus {{z,w) e C* X C* : P{z,w) = 0}. 

Gauge equivalent dimer weights give the same spectral curve. That is because after 
Gauge transformation, the determinant multiplies by a nonzero constant, thus not changing 
the locus of P{z, w). 

A formula for enlarging the fundamental domain is proved in [6l [T7]. Let Pn{z,w) be 
the characteristic polynomial of Gn, and P{z,w) be the characteristic polynomial of Gi, 
then 

Pn{z,w) = n n ^(^'^)- 

u"=z v"-=w 



3.2 Duality Transformation 

We call a geometric figure built with a certain number of bonds a closed polygon if at every 
lattice point, only an even number of bonds occurs. We associate to each square grid a 
Fisher graph by replacing each vertex with a gadet, as illustrated in the following figure 





ai 








b 




a2 



\b2 



R 



as 



D 



Figure 2: vertex of Ising model and corresponding dimer gadget 

If we draw a larger picture we will see that the Fisher graph is exactly the one obtained 
by replacing each vertex of the hexagonal lattice by a triangle. There is a one-to-one 
correspondence between closed polygon configurations on the square grid and dimer con- 
figurations on the Fisher graph. An edge is present in the closed polygon configuration of 
the square grid if and only if the corresponding edge is present in the dimer configuration 
of the Fisher graph, see Figure 2. 
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There are two ways to associate an Ising spin system on a square grid to the closed 
polygon configurations on the square grid. As a result, an Ising spin system on a square 
grid can be associated to the dimer system on a Fisher graph in two different ways. These 
two different ways lead to a duality transformation, which plays an important role in the 
identification of the critical temperature, as we will see. Note that these correspondences 
holds when there is no external magnetic field. Namely, the Ising model is defined by ([T]), 
given /i = 0, in order to introduce the correspondences between the Ising model and the 
dimer model. 



3.2.1 Correspondence 1: Measure-Preserving Correspondence 

An Ising model on a square grid is associated to closed polygon configurations on the dual 
square grid as follows: if two adjacent spins are of the same sign, the dual edge separating 
them is not present in the closed polygon configuration, otherwise the dual edge is present 
in the closed polygon configuration. An example of this correspondence is illustrated in 
Figure 3: 



Figure 3: Measure-Preserving Correspondence 



The dashed line are bonds of the Ising model, and the black line are edges of the dual 
grid. If two adjacent spins are of the same sign, the dual edge separating them is not present 
in the closed polygon configuration; otherwise the dual edge is present in the closed polygon 
configuration. This correspondence is two-to-one because negating the sign of all spins will 
end up with the same closed polygon configuration on the square grid. Combining with 
the one-to-one correspondence of closed polygon configurations on the square grid and the 
dimer configurations on Fisher graph, by assigning edge weights e~'^'^'= to the dual edges, 
we get a measure-preserving correspondence between the Ising spin system and the dimer 
system, namely, 

Pr{D{a)) = Pr{a) + Pr{-a) = 2Pr{a), 
where D{a) is the dimer configuration corresponding to the spin configuration a. 
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3.2.2 Correspondence 2: Correspondence Based on High- Temperature Ex- 
pansion 

Let ^ be a finite square grid, the high temperature expansion of the Ising model on Q is 

Zg,i = XI n exp(JeC7„C7^) 
<7 e=uveE{G) 

= X (cosh Je + au(Tv sinh Je) 

e=uveE{G) 

= ( JJ coshJe)^ JJ (1 + C7„(7„tanh Je) 

e=ui;eB(G) o- e=uveE{G) 

= ( JJ cosh Je) X Jl 2^*' tanh Je, 

e=uveE{G} CeS eeC 

where S is the set of all closed polygon configurations of Q. This way, up to a multiplicative 
constant, the Ising partition function is the same as the partition function of the closed 
polygon configuration of the same square grid. There is no one-to-one correspondence 
between configurations in this case, hence although in this case the partition function is 
invariant, it is not measure-preserving. Combining with the correspondence between the 
closed polygon configurations on the square grid and the dimer configurations on the Fisher 
graph, we obtain a correspondence between Ising models on square grid and dimer models 
on Fisher graph, by giving the corresponding edge of the Fisher graph tanh Je. 



3.3 Duality Transformation 

Let Gn be the quotient graph of the square grid on the plane, as defined on Page 4. Let 
G* be the dual graph of Gn- Define an Ising model on G„ with interactions {Je}eeE{Gn)- 
Assume the Ising model on Gn has partition function Zq^,!- Then Zq^,! can be written 
as, up to a constant multiple, following from the measure-preserving correspondence, we 
have 

^G„,7 = 2 n ^^P('^-) Yl lleM-'^Je)--='2 n exp(Je)ZF„,Ooo, 
eeE{G„) C''eS*oeeC* eeE{G„) 

where Sqq is the set of closed polygon configurations of G* , with an even number of occupied 

bonds crossed by both 7x- and jy. The sum is over all configurations in 5gQ. Similarly, 
we can define 5'oi('S'|o! •S'li) to be the set of closed polygon configurations of G*, with an 
even(odd,odd) number of occupied bonds crossed by jx, and an odd(even,odd) number 
of occupied bonds crossed by jy. Fn is the Fisher graph obtained from G* by replacing 
each vertex with a gadget, as described in Figure 2. Let Zp^^D be the partition function of 
dimer configurations on F„, with weights e~^'''= on dual edges of G„, and weight 1 on all 
the other edges. Then 

where 

^F„,De,, = XI n exp(-2Je). 

C*e5*^ eeC* 

For example, ^f„,_Doi the dimer partition function on Fn with an even number of occupied 
edges crossed by 7a; , and an odd number of occupied edges crossed by 7y. It also corresponds 



16 



to an Ising model which has the same configmation on the two boundaries parahel to ^y, 
and the opposite configurations on the two boundaries parallel to ^x- Similar results hold 
for ah the Zp^^Dg^ri G {0, 1}. 

On the other hand, if we consider the high temperature expansion of the Ising model 
on Gn, we have 



,/ = ( n ) E n 2"' tanh Je , 

e=™eB{G„) Ce5eeC 

where S is the set of all closed polygon configurations of G„. Let Fn be a Fisher graph 
embedded into an n x n torus, obtained from G„ by the correspondence in Figure 3, with 
weights tanh Je on edges of Gn, and weight 1 on all the other edges. In other words, the 
edge with weight tanh Je of Fn and the edge with weight e~^'^'= of Fn are dual edges. Then 
we have 

^G„,/ = 2 JJ exp(Je)ZF„,Doo = 2" n coshJeZ^^^^, 

ee£;(G„) ee£;(G„) 

where Zp ^ is the partition function of dimer configurations on Fn- 
Similarly, we can expand all the Zp^^^g^ as follows: 

ZF„,De,. = ^ n (l+exp(-2Je))Z^^^^((-ir,(-l)^). (15) 

Zp ^(—1, 1) is the dimer partition function of F„ with weights of edges dual to the edges 
crossed by multiplied by —1. Similarly for Zp ^(1,-1) and Zp ^(—1,-1). To see 
why (15) is true, let us consider, for example Zp^^u-^^, an odd number of present edges are 
crossed by ^x, and an even number of present edges are crossed by jy. The corresponding 
Ising model has boundary condition such that winding once along 'jx, the spins change 
sign, while winding once along 7^, the spins are invariant. This is the same as for a row of 
edges crossed by 7^, we change the coupling constant from Jg to — Je, and give the periodic 
boundary conditions. After the duality transformation, the edge weights tanh Je, if e is 
crossed by 7^ will be the opposite. 

Without loss of generality, assume n is even. Let K{z,w) be the Kasteleyn matrix of 
the toroidal graph, as defined on page 4. Given the orientation in Figure 1, we have 

II (l + exp(-2Je)) 



{[-PiKp^ (1, 1) + PiKp^ (1, -1) + PiKp^^ (-1, 1) + PiKp^ (-1, -1)] 
-hPiKp^^ (-1, 1) + PfT^^^ (-1, -1) + Pfi^#„ (1,1) + Pf^^„ (1, -1)] 
-l-PfKp^^ (1, -1) + PiKp^ (1,1) + PiKp^ (-1, -1) + PiKp^ (-1, 1)] 
-[-PfiC^,^ (-1, -1) + PfKp^ (-1, 1) + PfKp^^ (1, -1) + PiKp^ (1, 1)]} 

= n (l + exp(-2Je))Pf7^^„(l,l). (16) 

ee-B(G„) 

If we perform the same expansion for all the Pf JCi?„((— 1)^, (—1)'^), for 0, r G {0, 1}, we get 
Pfif^„((-l)^(-l)-) = -^ n (l+exp(-2Je))Pfi^^J(-l)^(-l)-), if{0,r}/(O,O). 
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Proposition 3.1. Consider a Fisher graph Fn embedded into an nx n torus. If an edge e 
is parallel to 7x,7y, give e weight We = e~^"''=. For all the other edges give weight 1. Define 
the duality transformation of Fn to be another Fisher graph Fn, embedded into an n x n 
torus. If an edge e is parallel to jx, ly, give e weight tanh Jg. For all the other edges, give 
weight 1. Let P{z,w) = denote the spectral curve. Then 

p^J(-l)^ i-m = ^ PpM-i)', {-m = o. 



4 Critical Temperature and Spectral Curve 

The purpose of this section is to identify the critical temperature of the periodic Ising 
model with the condition that the spectral curve as a single real node on T^. Using 
the correspondence based on the high-temperature expansion, we express the spin-spin 
correlation of the periodic Ising model as the determinant of a block Toeplitz matrix. The 
analysis of the holomorphic properties of the matrix leads to the result. 

Lemma 4.1. (Kelly and Sherman ll If ) Given a Hamiltanonian 

= — ^ Jac^a, 
AdA 

with ferromagnetic interactions, that is 

Ja>0, 

for all A C A, then 

^^ = {aBac)-M{ac)>0. 

The high temperature expansion of the partition function gives 

Z = 2 cosh Je tanh Jg, 

e C eeC 

where the sum is over all closed polygon of the square grid. Assume Tg = tanh Jg. We 
consider the spin located at (0,0), fioo and the spin located at (0,A^), ctqn. Then 

(O'OOO'OJV) = (<700O'0l)(<70lO'02)---(o"0,Af-l<70Ar)- 

Let 6/ denote the edge connecting (0, /) and (0, 1 + 1). Since 

o'o,zO'o,j+i(l + re,(To,i(To,/+i) = re,(l H aoio-o,i+i), 

we have 

N-l N-1 ^ ^ 

=z i[]2coshJe n 2^ 

e 1=0 o- /=0 e 

where Y[' means the terms corresponding to e/, < / < — 1 are omitted. This expression 
is of the form of a partition function for some Ising counting lattice with bonds on the 
straight line connecting sites (0,0) and (0, A^), replaced by We transform the square 
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grid Gn to the Fisher Graph, using the technique as described in Figure 3. Assume the 
edge weights of the Fisher graph are tanh Jg for all edges corresponding to edges of Gn, 
and all the other edges have weight 1. 

Let 7x(7y) denote a path winding horizontally(vertically) once on the torus, multiply 
all the edge weights crossed by ^xily) with z or ^, {w or ^), according to the orientation 
of edges. 

Let Kmn{z,w) denote the corresponding weighted adjacency matrix. If both m and n 
are even, we have 



(o'OOO'OAf) 



pfx.:„„(i, 1) + pfKU(-i, 1) + pfi^L.(i, -1) + pfi^;nn(-i, -1) 

PfiC^„(l, 1) + ViKmn{-l, 1) + Pfi^mn(l, -1) + Pfi^mn(-1, -1) ' 



where K'^^{z,w) denote the corresponding weighted adjacency matrix by changing the 
edge weights Tg; to . 

Theorem 4.2. For the Ising model on the infinite periodic square grid, we can express 
limjv^oo(o'ooCoAf)^ the determinant of a block Toeplitz matrix multiplied by a function 
of weights of all edges on the straight line connecting aoo and (Joat . 

Proof. We consider pf^'""(„'^) where u,v E {1,-1}. Assume 5 = K'^^{u,v) — Kmn{u,v), 
then 5 is given by 



(5(0, /; 0, 1 + 1) = -5^(0, / + 1; 0, 



R 




L 





U 





D 


^0 



R L 

J 







U D 



Te 



\ 





0/ 



If < Z < AT — 1 and zero otherwise. Here we are only interested in the entries labeled 

by R, L, U, D, because they are the only vertices adjacent to the edges with changed edge 
weights (from Tg; to — ) when we compute the spin-spin correlation. Thus if we define y 
as the 2N x 2N sub-matrix of 6 in the subspace where 6 does not vanish identically, and 
if we define Q as the 2A" x 2^" sub-matrix of -fC^n ™ the same subspace, we find 



(-l)'=^Pf (5 + Kmniu, v))Pi{K-l{u, V)) = (-l)'=^Pfj/Pf (y-^ + Q{u, v)), 



where ki,k2 are constants depending only on rn,n,N, and independent of u,v, and y is 
given by 



0,0 0,1 ... 0,Ar-l 0,1 0,2 ... 0,Ar 
R R R L L L 



0,0 
0,1 

0,A-1 
0,1 
0,2 

0,N 



R 
R 

R 
L 
L 

L 














-T 
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so that 



N-l 
/-n ^^l 



N(N-l) 
2 



Hence we have 



where s{u,v) is +1 or —1 depending on u,v. We choose a subsequence m',n' oo such 
that ^2 are always even, then 



7V-1 



hm (aooo-ojv)^ „ = IT ~ ^eJ^ 



m ,n -^oo 



1=0 



If we use K to denote the weighted adjacency matrix of the infinite graph, we can express 
Tjv as follows: 



1,1,K;0,0,K 



Q,N -1,R;Q,0,R 

K-^ 4- ^"1 

-^0,l,t;0,0,-R 

-1 °1 



0,0,R;0,Af — 1,R 

o,i,R;0,w— i,R 



0,l,-L;0,Af-l,R 



0,2,I.;0,Ar-l,R 



V 



0,JV,i:,;0,0,R 



^0,iV,I,;0,N-l,fl 



^0,0,BiO,l,i i_^2 
-1 1 
^0,1,R;0,1,L 



0,Ar-l,R;0,l,L 


ft— 1 

0,2,i:.;0,l,i 



^O.iV.LlO.liR 



■^0,0,R;0,iV,j:, 



0,1,R;0,N ,L 



0,N-1,R;0,N,L 

0,l,i;0,Af,Z. 

ft— 1 

0,2,I,;0,iV,I, 



It is well-known that Tjv is independent of u,v, since the entries of K ^ can be expressed 
as follows 



K ^ij,k,p;j',k',< 



1 

4^ 



2tt i-I-k 



JO 



^i<^l (fc-fc') g'i(/)2 O'-i') J^-^l (^g**^! ^ g«02 



(17) 



where K\\{z^w) is the weighted adjacency matrix of one fundamental domain embedded 
into a torus. And 



lim ((7oo(Toiv)^ = lim detTjv TT (1 - Tg^,)^. 



(18) 



0<KAr-l 



Let Zo denote the length of one period, assume AT is a multiple of Iq- We rearrange the 
rows and columns of T/v, in such a way that the determinant of T/v will not change after 
this manipulation, neither will the limit as A — t- oo. Then we group the matrix after row 
and column rearrangements, denoted still by Tjv, into ^ x ^ blocks, and each block is a 
2lo X 2Zo submatrix, such that the (j, /) block of Tj^ has a row range from (0, {j — l)lo, R) to 
{0,jlo — l, R), from (0, {j — l)lo + l-L) to (0, j7o, L), and a column range from (0, (/ — l)Zo, R) 
to (0,/Zo - 1,R), from (0, (/ - l)/o + l-L) to {0,flo,L). 
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If we define V'l a matrix- valued function on the unit circle, as follows 



1 

2^ 



2-K 



K-HC,e''^^)0,l,L;0,0,R+^: 



K-HC,e''^no,lo,L;0,o,R 

0,0,i?;0,l,L — i_^2 



K-HC,e''f'^)o,o,R;0,io- 



1,R 



0,lo-l,R;0,lo-l,R 



0,l,L;0,/o-l,R 



^-nC,e^*^)o,/o,L;0,Zo-i,R+T^ 



i^-nC,e''^^)o,Zo-l,i?;0,l,L 
i^-nC,e^*^)o,l,L;0,l,L 



^~nC,e^*^)o,0,i?;0,«o,L 
i^-nC,e^'^^)o,io-l,K;0,/o,L-i5|- 

i^-nC,e*'^^)o,l,L;0,/o,L 



\ 



i(f>2' 



0,lo,L;0,lo,L 



J 



then the block of "0 is the (j — j') Fourier coefficient of ip, given (17). Let T denote 

the infinite matrix as — )• cxd. Namely, 



(Tn) 



1 

2^ 



ICI=i < 



where the integral is evaluated entry by entry. Let T denote the infinite matrix as — )• oo, 
then T is a block Toeplitz matrix with symbol ^. □ 



Lemma 4.3. ({2^) LetT^ = {z, 



w z 



1, \w\ = 1} be the unit torus, either det K (z , w) 



has no zeros on T , or it has a single real node on T . That is, det K{z, w) = intersects 



at one of the points (±l,ibl), and both Q^^d 



ddet K{z,w) 
dw 



vanish at that same 



point. 



To study the asymptotic behavior of (o-Qoaoiy) as iV — >• oo, we first introduce the 
following lemmas 

Lemma 4.4. (WidomJW^) Let Tn{ip) be a block Toeplitz matrix with symbol ip, and n x n 
blocks in total. If 



)2 < oo 



fc= — oo 



k=—oc 



det i;{e'^) + 0, — Ao<e<2. arg det i^ie'^) = 0, 
Ztt 

here ||^fc|| denotes the Hilbert- Schmidt norm of the matrixipk, cf^c^ 27r^o<e<27r arg det 7/^(6*^) 
is the index of the point with respect to the curve {det ^/'(e*^) : < < 27r}, namely, 



1 1 

— Ao<e<27rargdetV'(e*^) = — 



1 



adetV(C) 



ICI=i 



det^(C) dC 



dC, 



then 



lim ^5^^=ii;[^], 
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with 

Gi'iP] = exp{^^'"logdetV'(e*')d^}, 

E[tp] = det T['ilj]T[tp-^]. 

where the last det refers to the determinant defined for operators on Hilbert space differing 
from the identity by an operator of trace class. 

Lemma 4.5. For ferromagnetic, periodic Ising model with nearest neighbor pair interac- 
tions on a square grid, if det K{z,w) has no zeros on T^, then 



lo 



1=1 

for any 6 G [0, 2tt). 



Proof Let us consider an m x n torus, where m and n are the numbers of periods along 
each direction. The previous technique to compute the spin-spin correlations implies that 

7' 



l<«</o 



mn 



where is the partition function of dimcr configurations on an m x n Fisher graph with 
weights Te((l <l<n) changing to — . Since the graph is embedded into an m x n torus, 

(Too and uon are actually the same spin, therefore we have (aoocron) = 1, that is 



n 

l<l<lo 



'ei^mn ~ ^mn- 



Let K'^J^Kmn) be the corresponding Kasteleyn matrix with respect to Z'^J^Zmn)- We 
are going to construct a correspondence between monomials in r^^ det ^(1, tu) and 
det iir^^i(— 1, w\ where w is any number on the unit circle. Each term in the expansion of 
the determinant corresponds to an oriented loop configuration in the m x 1 torus, where 
each vertex is adjacent to exactly two edges in the configuration. The correspondence is 
constructed by changing configurations on the first row of the graph, while keeping the 
configurations on all the other rows unchanged. Consider the first row. For each pair of 
nearest components, which are parts of non-doubled edge loops crossing the first row, there 
exists a pair of nearest triangles, one on each component. We change their paths through 
the pair of nearest triangles. This way, in the new configuration, all the ei occupied once 
originally will still be occupied once, with weight ^ changed to Tg, . Any ei unoccupied will 

be a doubled edge, with weight 1 changed to Tg^ any e/ occupied twice will be unoccupied, 
with weight changed to 1. Hence the corresponding terms in ni<KZo ^e, det w) 
; r 1 1 1 det i^m,i( — 1, w) have the same absolute value. They arc equal if the signs are also 
the same. For each component, which is part of a non-doubled-edge loop crossing the first 
row of the torus, if we change the path through a pair of triangles (corresponding to two 
nearest components on both sides), the sign will change by a factor of -1. This sign change 
will cancel if we change the sign of the z-edge passed by the loop simultaneously. Hence 
we have 



n 



^li deti^;„,i(l,u;) = detiC„,i(-l,u;) 
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An example of such a path change is illustrated in Figure 4. Similarly we have 



II TldetK'^^i{-l,w) = detKm,iil,w). 

l<l<lo 



Therefore, we have 

lo 

llld- 

1=1 



(n|°=ire^,)detK;,,i(l,e^^) 
detK^,i(l,e^^) 

detir^i(-l,e^'') 



detV'(e" 



lim 

m— >-oo 



lim 



m^oo det Kmi{l,e^^) 

Let vi, - ■ ■ , V2nio be vertices which are endpoints of those edges changing weights from 



to — . Let us consider 



detK'^^,{l,e^<^) 
m-^oo del Km,i{l, e*^) 



lim 



,Too 2- detK„i(l,e^^) 



2n 



Ul,--- ,U2nlQ,Uj~Vj j=l 

where K'^^{l,e'^)Ec{ui,- ,u2ni^} is the submatrix of K'^i{l,e^^) deleting rows vi,--- ,V2nio 
and columns ui, - ■ ■ ,U2nio- '^^^^ submatrix is also a submatrix of XmiCl, e'^), because all 
the edges with changed weights have been deleted. K!^i{l, e^^)yj^uj are edge weights, which 
arc independent of the first variable - weights of edges crossed by 7,,.. As m — t- 00, the limit 
of entries of Kl^\{l,e''^) is independent of the first variable too, since liiUm^oo K^\{z, e'^^) 
can always be expressed as the same complex integral no matter where z is located on the 
unit circle. As a result, 



detK^i(l,e*^) m-)-oo detifTOi(— l,e*^)' 



Hence we have 



lim 

m->oo 



detKmi{l,ei^) 



.. detKmi{-l,e'') 

lim ^— 

m^oo detiCmi(l) e*'') 

detjir^i(l,e^^) 
m-*-oo detif^i(— l,e'^) 

nti<deti^;,i(-l,e'^) 
m->-oo detKTOi(— l,e*^) 



Therefore 



[n(l-rfj2detV(e^ 



= 1. 



(19) 



When 9 = 0, the quantity is the quotient of two skew-symmetric matrix, hence it is positive 



and equal to 1. The lemma follows from the continuity of [H^Lill ~ '^tej^ ^et ipie 



□ 
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Lemma 4.6. (Gohberg,Krein^lO]) Let denote the set of all completely continuous op- 
erator A for which 

oo 

^Sj{A) < oo, 
i=i 

where sj are singular values of A. Let A(fj,) be an operator-function with values in and 
holomorphic in some region. Then the determinant det(/ — A{fi)) is holomorphic in the 
same region. 

Lemma 4.7. If we allow complex edge weights, then the spin correlation 

lim (o-ooo-oAr)p 



under toroidal boundary conditions is analytic with respect to edge weights except for the 
edge weights tanh Jjj- such that the spectral curve has a real node on T^. 

Proof. Let Hlip] denote the semi-infinite Hankel matrices with kernel ■i/', namely, 

-ffiV'] = (V'i+i+i), 0<i,i<oo, 

then 

T[,lj]T[^~^\=I-H[m[i^~\ 

where 

Since "0 is an even-order, anti-Hermitian matrix on the unit circle, its determinant will 
always be real, hence we have 

Ao<0<27rargdetV'(e*^) = 0. 

For any edge weights which do not lie on the critical surface, the condition of Theorem 4.4, 
and 4.5 will satisfy(Widom^28j). Then we have 

lim (aoocTo7v)p = lim TT (1 - r^f detT,v[^] = det(/ - H[^]H[^-^]). 

0<i<Af-l 

Since the entries of ii\\\)\H\\\)~^\ will be analytic everywhere except on the spectral curve, 
so is the spin correlation limAr_5.oo(o'ooO"oAf)p- D 

Lemma 4.8. The weak limit of the Boltzmann measures Vn on Gn, the quotient graph 
defined on page 4, is a translation invariant Gibbs measure V. The probability of occurrence 
of a subset of edges ei = uivi, Cm = UmVm of G in a dimer configuration of G chose 
with respect to Gibbs measure V , is 



P(ei,...,e„) = J]z/;,,,jPf(i^ 



'ei,...,em I ) 



i=l 



where K ^ is the inverse of the infinite Kasteleyn matrix, defined as in and {K ^)ei,...,e, 
is the sub-matrix of , whose lines and columns are defined by the vertices of the edges 

6l ; • • • ; Cm . 
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Proof. The dimer partition function Z„ of the graph G„, for n even, is 

Let E = {ei = uivi, ■■.,em = UmVm} be a subset of edges of G„, then the probabihty 
P„(ei, ...,em) of these edges occurring in a dimer configuration of G„ chosen with respect 
to the Boltzmann measure Vn, for n even, is 

nm 



where i?'^ = V{Gn) \ {uivi, ...,UmVm} and {K^)e<^ is the sub-matrix of whose hues 
and columns are indexed by E"^. Let P(z,w) = det K{z,w). By lemma 4.3, either the 
spectral curve P{z,w) = does not intersect T^, or it has a single real node on T^. 
If P{z,w) = does not intersect T^, 



n-ln-l , , ^ rtr^, i^(2j + e) t7r(2fc + r) 

i(2j + 0)7r{x-x') i(2k + T)7v(y-y') C Q f [K { 6 ri ,6 " )\y y' 

e " e 



^-^n ){x,y,v),(x',y',v') „2 5Z 5Z o - i7r(2j + fl) i7r(2fc + T) 

j=Ofc=o P{e " ,e n j 

The right side is a Riemann sum in non-critical case, it converges to an integral for any 
0,T £ {-1, 1}, that is, 

J_ ff ^x~x'^v-v' Cof[K{z, w)]yy dz dw 
4vr2 J Jj2 P{z, w) iz iw 

We also have 

pf«v=(-i)'^^"^^^^S^, 

n 

where k : N — t- {0, 1} is a function. We can choose a subsequence {n'} such that n{n') is a 
constant kq. Then 

P(ei, ...,em) = lira Vn'{ei,...,em) 

T-rm 

n'— >oo -^^n 

+Pf(i^l?)^lpfK^? + Pf(K^l)^lpfi^ll]| 
m 

i=l 

Note that the convergence actually holds for all n, because on the complement subsequence 
N\{n'}, K is also a constant given that ac is a two- valued function. The sign does not matter 
since we take the absolute value in the end. 

If P{z, w) = has a single real node on T^, see [Ij for a proof. The weak convergence 
of Vn for n G N follows from the fact that there is a unique translation invariant Gibbs 
measure, and the weak limit of Vn on any subsequence is translation invariant. □ 

Lemma 4.9. Let V be the weak limit of Vn, Boltzmann measures defined for graph Gn 
with periodic boundary conditions. P{z,w) = does not intersect T^, if and only if the 
edge- edge correlation 

V{eike2) -V{ei)V{e2) 
converges to exponentially fast as |ei — 62] oo. 
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Proof. By Lemma 4.8, 



V{eike2) -V{ei)V{e2) 



ViK- 



Ui Vi U2 V2 
Ui Vi U2 V2 



uivi Wu2V2 



U2,V2 



K-^ I 

-"■Ml,n2-"-i>i,i'2 I 



(22) 

■ \K'^ I 
(23) 



If P{z,w) = has no zeros on T^, (21) and (23) imply that as |ei — 62] 
converges to exponentiaUy. If P{z,w) has a real node on T^, then 

Cof{K{z,'w))yy 



00, P(ei&e2) 



P{z,w) 

is not an analytic function in edge weights, hence its Fourier coefficients does not decay 
exponentially fast. To see that, if the Fourier coefficients decays exponentially, then the 
function is the limit of a locally uniform convergent power series, hence the function has 
to be analytic. □ 

Lemma 4.10. If P{z,w) 
det ( 



has a single real node on T^, and if 



902 amp 



(e,<i>)={m 
then the edge-edge correlation 

V{eihe2)-V{ei)V{e2) 
converges to quadratically as |ei — 62] — s- 00 

Proof. Without loss of generality, assume the real node is (1, 1). Then 

= P(l,l) = (Pfi^(l,l))2, 

where PfX(l, 1) is linear with respect to each edge weight. Fix an edge e = uv, and assume 
the weight of edge e is We, then 

P{1,1) =C2wl + CiWe + Co, 

where C2 is the signed weighted sum of loop configurations in which e is a doubled edge. Ci 
consists of two parts: Ci = Di + 'D2. Vi is the signed weighted sum of loop configurations 
in which e is part of a loop oriented from u to v; T>2 is the signed weighted sum of loop 
configurations in which e is part of a loop oriented from v to u. When z = w = 1.,'Di = 'D2. 
Co is the signed weighted sum of loop configurations in which e does not appear at all. The 
determinant of the submatrix of iC(l, 1) by removing the row indexed by u and the column 
indexed by v is exactly C2We + T>i. 
Since P(l, 1) is a perfect square 



C0C2, 
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as a result, 



P(hl) 



{C2We + Vi)We + {ViWe + 



C2' 



1 



{C2We + Vi] 



P(l, 1) = implies C2We + T^i =0, hence the determinant of a submatrix of K{\, 1) by 
removing the row indexed by u and the column indexed by v is 0, if uv is an edge. Now, 
let us consider AdjK{l, 1), i.e. the adjugate matrix of K{1, 1) whose entry is (— 
multiplying the determinant of the submatrix of K(l, 1) by removing the jth row and the 
ith column. Since deti^(l, 1) = 0, AdjK{l, 1) is a matrix of rank at most 1. As explained 
above, the entries of AdjK{l, 1) corresponding to an edge is 0. This implies that AdjK{l, 1) 
has entries on each row and each column. Therefore AdjK{l, 1) is a zero matrix, and 
any minor of 1) by removing a row and a column is 0. 
In a neighborhood of {9, (j)) = (0, 0) P has the expansion 



ar + 



where a, 6, c G R are real numbers and • • • denotes terms of order at least 3. Then near 
a node, a point on P = satisfies either 9 = Xcf) + 0{(j)'^) 01 9 = X(j) + 0{4P'), where A, A, 
which are necessarily non-real, are the roots of ax^ + 6x + c = 0, since P{z, w) > 0, for any 
{z,w) G \ {(1, 1)}([24J), and the assumption that the Hessian matrix at (^,1;^') = (0,0) 
is invertible. This implies 



4ac < 0, 



d^P 



89"^ 



>0, 



d^p 



(e,</.)=(o,o) 



> 0. 



(e,0)=(o,o) 



Since P{z, w) 
obtain 



Aet K{z,w), we apply the differentiation rule for the determinant and 



dP 
~d9 



^ det 



where is a matrix which are the same as K except the row indexed by v. For each entry 
on the row indexed by we take the derivative with respect to 9. The sum is over all the 
vertices of the quotient graph. Then we expand det with respect to the row indexed 
by V. Obviously at most one entry on that row of Ky is nonzero. Hence det Ky = or 
daiKy = ±ie^''^{AdjK) y^y), where vw is an ^;-edge or ^-edge. Then we take the second 
order derivative at 6* = 0, (/> = 0, 



d det Ky 
89 



e^''{AdjKil,l))y,y,±e 



±ie d{AdjK)y. 



d9 



(e,0)=(o,o) 



The first term is 0. However, since 

d{AdjK)^^n, 



/ 0, there exists at least one z-edge vv' 



such that 



(e,0)=(o,o) 



(e,<^)=(o,o) 

7^ 0. When [9, 4>) lies in a neighborhood of (0, 0), we have 



p9 + . 



P{e'' 



a9^ + 
1 

A- A 
1 



+ cq 
pX — 



+ 0(1) 



x<p) 

pX 



+ 



a{9 - A(/>) 



+ 0(1) 



+ 



q-pX 



A — A ia{z — 1) — aX{w — 1) a{z — 1) — aX{w — 1) 



+ 0(1), 



28 



where z = e^^ ,w = e^"^ . It is proved in [T7] that if 

R{z,w) = a{z - zq) + j3{w - w^) + 0{\z - zqP + \w - wq^'), 

then we have the following asymptotic formula for the Fourier coefficients of R^^ 

1 I w^z^ dz dw "^^o-^o + 0( ^ 

{27ri)'^ Jj2 R{z,w) z w 27ri{xazo — yl3wo) + 

Therefore we have 

1 f w^zyCof[K{z,w)]yy dz dw 

+ 0( ' • 



{2TTiy Jj2 P{z,w) z w 

1 r Q — _|_ — Q 



A- A 



2TTi{xa + yaX) 2TTi{xa + yaX) 



rj.2, _|_ yl 



Hence K^^, decays linearly as |u — v'| — >• oo, as a result, the edge-edge correlation decays 
quadratically as — f'| — oo. □ 

Lemma 4.11. Let Vi be a Gibbs measure of the ferromagnetic Ising model on the periodic 
square grid, V be the induced Gibbs measure of the dimer model by the Fisher correspon- 
dence in Figure 2. If under Vj the spin-spin correlation decays exponentially, then under 
V , the edge-edge correlation decays exponentially. 

Proof. Let pi = be the lattice gas variable. Let A, B be two finite subsets of the 
vertices of the square grid. Define 

PA = Wpi 

i&A 

Clearly Sa^ Pa and Sa — PA ai"e nonnegative, monotone functions of the and thus also 
of the {(Tj}. According to the F.K.G inequality, we have 

{pb{Sa -Pa))> {pb){Sa - pa) > 0, (24) 

and 

{Sa{Sb - Pb)) > {Sa){Sb - Pb) > 0. (25) 



Combining (24) and ( |25[ ), we have 

< {paPb) - {pa){pb) < (SaPb) - {Sa){pb) < (SaSb) - {Sa){Sb) 

ieAjeB ieAjeB 

Let 61(62) be an edge with endpoints ui, vi{u2, ^2)- According to the Fisher correspondence 
illustrated in Figure 2, we have 

7'(6l&62) -V{ei)V{e2) = {pu^l- Pv^)pU2{'^- PV2)) - {puA'^- Pv^)){pU2{^- PV2)) 

+ {p^^{l- Pu^)py^{l- PU2)) - {pvA^- Pu^)){Pv2{'^- PU2)) 
+ (P«i(l - Pvi)Pv2{'^ - PU2)) - (/5«i(l - Pv^)){Pv2{^ - PU2)) 
+ (Pt>i(l - P«i)/5«2(1 - ^^-2)) - {PvA'^ - Pmi))(P«2(1 - PV2))^ 
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and 



\{Pm{l - Pt,i)p„2(l - Pvi)) - (P«i(l - P»;i))(P«2(l - PV2))\ 

= \{{PU1PU2) - (P«l)(P«2» - {{PuiPuiPvi) - {Pui){Pu2Pv2)) 

-Hp ) - {PuiPvi){Pu2)) + {{PuiPvi Pu2Pv2 ) - {pUiPvi){pU2pV2))\ 

< M{PuiPu2) - {Pui){Pu2)\ +^\{P ) - {Pui){pV2)\ 

+ '^\{Pu2Pvi) - {Pu2){Pvi) \ + \ {P V2/\- 

The other terms of V{ei&:e2) — 'P{ei)V{e2) can be bounded in a similar way. Hence if the 
spin-spin correlation decays exponentially, the edge-edge correlation decays exponentially 
after the Fisher correspondence described in Figure 2. □ 

Lemma 4.12. Given interactions {Jk}k=i' when j3 increases from to +oo, j3c corresponds 
to the smallest (3 that the spectral curve for dimer configurations with weights tanh/3Jfc has 
a real node on T^. 

Proof. Consider the weak limit of Boltzmann measures of the Ising model under periodic 
boundary conditions, 

{aiaj)p = {aiaj)p - {cri)p{aj)p 

= ^{{piPj)p - {Pi)p{pj)p) = MptPj)p - 1- 

Since pipj is an increasing function with respect to each single variable a, {(7iaj)p is an 
increasing function with respect to /3. In particular lim„_5.oo(o"ooO'On)p is increasing in /3. 
There is a unique critical temperature /3c,p, for any given interactions {Ji,j}i,j, such that 

lim {aooaon)p = if /? < /3c,p, 

n— >oo 

hm {aoocron)p > if /3 > (3c,p. 

n— >co 

Hence at the edge weights tctnh. j^c^pJ ^ lii^n-^^cxo (<7ooC'On)p fails to be an analytic function of 
/3(complex value of (3 are also allowed). By Lemma 4.7, tanh/3c,pJ is a subset of the edge 
weights such that P{z,w) = has a real node on T^. By Theorem 1.1, the uniqueness of 
the translation invariant Gibbs measure implies that 

lim (crooo-o„)+ = lim {croocron)p- 

Therefore for /3 < /3c,p 

< lim {aooaon)+ - {(^00)+ < lim (fioocron) = 0, 

n— >oo n— ^00 

where the first inequality follows from Lemma 4.1. Hence for f3 < f3c^p, the spontaneous 
magnetization 

Ai, the quotient graph consisting of one period. Therefore f3c,p < /3c- Pc, as defined before, 
is the reciprocal of the lowest temperature such that m* = 0. 

It is proved in [21 [1] that for /3 < /3c, the spin-spin correlation decays exponentially, 
therefore after the Fisher correspondence, the edge-edge correlation on a graph with weights 
Q-'^PJ decays exponentially. We derive that /3c < /3c,p, because at j3c,p, P{z, w) = has 
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a real node on T^, the edge-edge correlation does not decay exponentially, see Lemma 
4.9. Combining with the previous argument, we have /3c,p = Pc- When /3 < /3c, the 
corresponding spectral curve with weights e~^^'^ does not intersect T^. By proposition 
3.1, the spectral curve of dimer configurations with weights tanh/3J does not intersect 
either, if /? < /3c- This completes the proof. □ 

Lemma 4.13. Assume the Ising model has period n x n, n is even. If the weights of all 
edges parallel to ")x)ly in the open interval (0, 1), and all the other edges have weight 1. 
Given an orientation to the Fisher graph as illustrated in Figure 1, the only possible real 
node of the spectral curve P{z, w) = on is (1, 1). 

Proof. It suffices to prove that under the assumption of the lemma, Pfi^(l, —1) ^ 0, 
Pfi^( — 1, 1) 7^ and Pfi^(— 1, —1) 7^ 0. Given an orientation as in Figure 1, 

PfK(-l, -1) = Zd,oo + Zd,oi + Zd,io - Zd,ii, (26) 

where Z^^i is the partition function of dimer configurations with an even number of occu- 
pied edges crossed by 'yx and an odd number of occupied edges crossed by ^jy. Zufio-, Zd,io, Zd,ii 
are defined similarly. The expression of P(— 1, — 1) follows from the fact that if we reverse 
the orientation of all the edges crossed by 7^ and jy, the resulting orientation is an crossing 
orientation, as defined on Page 4. Similarly, we have 

PfK(l,-l) = ZDfio + ZDfli- Zd,io + Zdm, (27) 
PfK(-l,l) = Zd,oo-Zd,oi + Zd,io + Zdm- (28) 

Zd,oi also corresponds to an Ising partition function with the opposite spin configurations 
along the two boundaries parallel to jy, and the same spin configuration along the two 
boundaries parallel to ^x- After the duality transformation, Z^^qi becomes, up to a constant 
multiple, the partition function of dimer configurations on the dual Fisher graph F, with 



weights of edges crossed by 7x multiplied by —1, see ( 15 ). Since the duality transformation 
does not change the partition function, up to a multiplicative constant, and after the duality 
transformation Zj^^q, Z^^i, Z^^iq, Z^^h are partition function of the dimer model on the 
same graph with the same edge weights except on the boundary edges, and Zd^q is the 
only one of them such that all the edge weights are positive, as a result, Zd^q is the biggest 
of {Z£)^gT-}e.Te{o,i}^ because all the others are partition function of dimer configurations of 



the same graph with negative weights. As a result of (26), (27) and (28), and the fact that 
each one of {ZD,9T}e,T£{o,i} is positive as the partition function of dimer model on a graph 
with only positive edge weights, 

PfE:((-l)^(-l)") >0, if (0,r)/(O,O). 

□ 

Lemma 4.14. Assume n is even. For any given { Jfej^^a, < < 00, the curve 

■j{t) = (tanhtJi, • • • ,tanhtJ2„2) < t < +00 

in the edge weight space(a 2n^ -dimensional vector space with coordinates given by edge 
weights) intersects P{1, 1) = at a unique point j{to), < to < +00. 
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Proof. First of all, we claim that 

7(t)nmi,i) = o}/0. 

Note that P(l,l) = (PfK(l,l))2. Given edge weights 7(0) = (0, • • • ,0), Pfi^(l,l) = 1. 
Consider the edge weights limi_>_|_oo 7(t) = (!,••• , 1). After the duality transformation, all 
the edges parallel to 7x, have weight 0. Since ^d,01) -^d.iO) -^d.ii corresponds to dimer 
configurations by negating the weights z—{w—,z— and w—) edges, and all such edges have 
weight 0, we obtain 



for weights (!,••• , 1), see (15). Hence limt_!>+oo Pfi^(l, 1) < 0. There exists < to < c«, 
such that 7(to) G {^"(1, 1) = 0}. 

Now we prove the uniqueness of the intersection. Since -P(l, 1) = is invariant under 
the duality transformation, i.e, any edge weights lie on the surface -P(l, 1) = if and only 
if the images after duality transformation lie on the surface -P(l, 1) = 0, it suffices to prove 
that the image of j{t) 

7*(t) = (exp(-2tJi),--- ,exp(-2tJ2„2)) 

intersects P{1, 1) = at a unique point. Assume f{t) is the function obtained by plugging 
edge weights (exp(— 2t Ji), • • • , exp(— 2t J2„2)) into Pfi^(l, 1), and g{t) = 2 n^^i exp(t Jfc)/(t), 
it suffices to prove that there is a unique to in (0,+oo) such that 5(^0) = 0. Given inter- 
actions (tJi, • • • ,tJ2„2) on edges of the quotient graph of Gn, let Z^'^ j be the partition 
function of the Ising model on Gn- Let Zq^ li^c^n /> i) partition function 

of an Ising model, with the interactions on edges crossed by ^x{ly, Ix and 7^^) multiplied 
by —1. Then we have 

9^^) = ^GnJ ~ ^Gn,I ~ ^Gn,I " ^Gn,7' 

e=uv^E{Gn) e=uveE{G„) 

~^Gn,I X] WuCTv) 10 - Z}.\ J ^ Je^icTuCTv) 11, 

e=uv£E{Gn) e=uv£E(G„) 

where {^e'^}eeE(G„) are corresponding interactions for Zq^ j and ( )eT- are the expected 
values with respect to Boltzmann measure defined by {Je^}eeE{Gn)- ^0 be the subset 
of E{Gn) consisting of edges intersecting neither nor 'jy, and Ex{Ey) be subsets of E{Gn) 
consisting of edges crossed by '^x{ly)- 

9'{t) = ^ 'Je^{Z^„j{(Tu(yv)()0- Zq^^j {a u(J v)oi- Z^^ j {a uCfv) w- ZQ^^ i{auCF 

e=uv&Eo 

+ ^ Je^i^G'n,lWuCrv)oo — ZQ^^ j{auCrv) 01 + Zq^^ j {a uCTy) 10 + Zq:^ j {a uCTv) 11) 

e=uv(iEx 

+ X/ '^e°(^Gi,/(^«'^f)00 + Zq^^ I {a uCfv) 01 - Zg^„jWu(^v)iO + ■^G„,/('^«^f 
e=uv&Ey 

We claim that for {9, r) 7^ (0, 0), 

\{cruO'v)eT\ < {o'uO'v)oo- (29) 
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First of all, 

{cruO'v)eT < (o"«cr^)oo, (30) 



by F.K.G inequality. The inequality (30) is strict because the graph is finite. Moreover, 



Gn,I G„,I o-,(t' ee-B(G„) 

= ^00 ^^er ^(l + t^^.)^k^^.exp( ^ (jf + jOO))^^^^]^ 

G„,/ Gn,/ t o- ee-E;(G„) 

where we have introduced the Ising variables a'^ = ±1 an tj = a'^ai. Since Jg° + jf"^ > 0, 
the ferromagnetic condition and finiteness of G„ imply 

{(yu<yv)eT + (o"«'7i))oo > 0. (31) 



(l29j) follows from (jSOj) and (jSlj). Hence 

g{t) > (Zq^ j - Z^\ j - Z(9^ j - Zq^ j) ^ Je°{au<yv)oo- 

e=uveE{Gn) 



Since J°° > 0, (o-«fT„)oo > 0. 

Z'^'^ — Z^^ — Z^'^ — Z^^ = Pf-^(1; 1) 

Gn,I G„,I G„,I G„,I r,i-r2n2 u t \ ' 

2nfc=iexp(iJfc) 

We have if Pfi^(l, 1) > 0, g'{t) > 0, g{t) is strictly increasing as t increases. Since the 
duality transformation changes the sign of Pfi^(l,l), we have g{0) < 0, see (16), and 
liuit-^oo git) > 0. When t increases from to oo, originally g{t) is negative, and there 
g{t) = is possible for some t. Consider the first such t that g{t) = 0. From then on, g{t) 
will always be positive since g'{t) > 0. As a result, there is a unique to, such that g{to) = 0. 
This proves the lemma. □ 

Proof of Theorem 1.1 Theorem 1.1 follows directly from Lemmas 4.11, 4.12 and 4.13. 

Example 4.15. (1x2 Ising model) Consider an Ising model on a square grid whose 
interactions have period 1x2, as illustrated in Figure 5. The critical reciprocal temperature 
/3c is the solution of the following equation: 

P{z = l,w = -1) = detK(l,-l) 

= [62 + C162C2 - bib2 + C1C26162 - 1 + C2C1 + 61 + ci6iC2]^ = 0, 

where 



bi = tanh(/3cJ4), ^2 = tanh(/3cJ3), 
ci = tanh(/3cJi), C2 = tanh(/3cJ2)- 

In this case, w = —1, because the size of the period is oddxeven. 
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Figure 5: 1 x 2 Ising model and Fisher graph 

5 Simple Proof for the Exponential Decay of Spin Correla- 
tions with Symmetric Interactions 

The goal of this section is to prove the exponential decay of spin-spin correlations for 
symmetric interactions above the critical temperature. Namely 



-a\i-j\ 



where a > is a constant. The symmetric interactions are those who are invariant under 
the reflexion by the horizontal axis y = and the vertical axis x = 0, see Figure 6. Lemma 
5.1 proves the exponential decay along the lattice direction, and Lemma 5.2 proves that the 
spin-spin correlations along arbitrary directions are dominated by those along the lattice 
direction. Hence the result follows. 



Lemma 5.1. Under the assumption that det K{z,w) ^ 

(fOOCTOAf)^ 

where < K,a < oo are positive constants. 



lim (ctooo-oat)^ < Ke 

N-^oo 



Proof. According to (18), it suffices to consider the asymptotic behavior of Y\ 



0<l<N-iy 



TgJ^ det T/v[V']) as N large. This is the determinant of a block Toeplitz matrix, obtained 
from Tjv by multiplying the jth row (1 < j < A^) by (1 — Tg_,_^)^. The determinant of its 
symbol is 1, according to Lemma 4.5, since Te^ changes periodically. In the proof of this 
lemma, we denote the new Toeplitz matrix by the same notation T/v, and the new symbol 
by the same notation ■0, which is the original ip multiplied by a constant for each row. 
Assume 



lim Tn[ip], 



lim Tr,. 
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Define A,- P] Kr the Banach algebra under the norm 

oo oo 



fc=— oo fc=— oo 

where </) is an r x r matrix- valued function. || • || on the right side is any norm for r x r 
matrices, and 0^ is the A;th Fourier coefficient of (j). First consider the case that the 
Toeplitz operator T[ip] is invertible, which is equivalent to the matrix- valued function ijj 
has a factorization 

where ^/J± are invertible in Ar f] , and ^+ and ip- have Fourier coefficients that vanish 
for negative, resp. positive indices. The Hankel operator Hlip] is Hilbert-Schmidt, the 
product of any two such is of trace class. 

H[^] = (V'i+j+i), 0<i,j<oo, (32) 
r„[^(A]-r„[V']T„[</.] = PnH[i>]H[4>]Pn + QnH[i^]H[4>]Qn, (33) 

where P„ and Qn are defined by 

Pn(/0,/l,-) = (/0,-,/n,0,...,0), 
Qn(/0,/l,-) = (/n,...,/0,0,...,0). 



If we replace (f) by in (33), since H[%Ij_'^] = 0, we have 
Multiply the above equation by T„['(/'I[;^], we have 



Use -0+1 '0+"'^ to substitute in (33), we have 

/„ - T„[V.+]T„[0;i] = PnH[i;+]H[i^-^]Pn + QnH[iP+]H[ip-^]Qn 
Since H[^jJ^^] = and H[tjj^] = 0, we have 

Tn[^P+]Tn[i^+^] = In. 

Therefore 

T„[V']r„[0Zl]rn[0;'] = In - PnH[^|;]H[i^Z^]PnTn[^^'] 

= Pn{I - H[xl,]H[^Z^]PnT[i,Z,^])Pn. 

Since T„['0^^] are block triangular matrices, one sees that the left side has determinant 
exactly 

Dn[mr-'r'G[^-^'r^ = 

where -Dn[V'] is the determinant of r„[^], and G is defined as in Lemma 4.4. Moreover, 
explicit computations show that 

H[4,]H[rJ]T[r+^]=H[il;]H[^-\ 
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Let II • 111 denote the norm of operators in S^i(see Lemma 4.6, and the appendix). The 
uniform norm | • | of an operator A is the number 



sup |^^|. 
|<^|=i 



We have 



\\PnH[^j]H[ir']Pr,T[^-']Pn - H[^/j]H[^^'^]\\l 

< \\PnH[i>]H[i;Z']PnT[i>+^]Pn - PnH[iP]H[ipZ^]T[,p^^]P4i 
+\\PnH['il;]H[^-^]Pn - H[^P]H[4>-^]\\i 

< Ke-"", 

where < K, a < cx) arc positive constants. To sec the last bound, note tliat H[il)Z^]Pn — 
H[^Z^] has the form ^ || * ) ' ^'^^^^ nonzero entries correspond to columns of H[ijjZ^] 

with indices greater than n, which are exponentially small by definition, similarly for 

PnHlijj] - H[i}] .Define 

An = PnH[i;]H[4^Z^]PnT[i;-^]Pn, 
A = H[^]H[i;-\ 



then 

Define 



Pn-^||i < iiTe""". 

DaM = det{I-nAn), 
Da{i^) = dei{I-nA), 
A{fi) = A{I-fiA)-\ 

Anin) = An{I-nAn)~^. 

We shall call the complex number /i an F-regular point of the operator A, if the operator 
/ — fiA has an inverse. Let F be a simple rectifiablc contour, consisting of F-rcgular points 
of the operator A, which encloses the point = and the point = 1. By virtue of the 
maximum modulus principle, if we can prove that for any m„ — A||i < Ke~°'"', 

max|Z?A(M) - DaM\ < Ci^e-"", 

then 

\DA{l)-DA„{l)\<CKe-''^. 

Let us denote by L a simple rectifiable curve consisting of F-rcgular points of the operator 
A, which connects the point /x = with some point of the contour F. Let F^(/x G F(JL) 
denote the shortest path along the curves F and L which connects the point /x = to the 
point /i. From the equality 

(/ - Mn)"' = {1- - KAn - A){I - flA)-']-\ 

it follows that when the condition 

\\An-A\\^< min [H|(7-M)-i|]-i 
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is fulfilled, all points fi of the curves F P| L are F-regular points of the operator An and 

max \(I - LiAny-^l < C max \(I - llA)~-^\, 

where C is a constant depending only on the curves F and L and the operator A. Since 
Ail,) - AM = (/ - fiA)-\A - An){I - f,An)-\ 

we have 

\\A{fi) - An{fi)\\l <\{I- - ^nlllKl - Mn)-'|. 

for any fiQ G F |J L. Then 



tr{A{^l) - Anifi))dfi 



< 



\\A{i^)-Anii^)\\idfi<C\\A-An\ 



where C is another constant depending only on the curves F and L and the operator A, 
and tr{K) is the sum of all eigenvalues of operator Er(for details, see the appendix). F^q 
is the shortest path along the curves F and L connecting the points and /io- Then 



1 



Da{^^q) 



and 



log 



DaM DAils) 



r(A) 



tr{A{ij)-An{ii)). 



where r{K) is the dimension of operator K, namely, the dimension of the closure of the 
range of K. Hence 

\DA{^l) - DaM\ = I^a(/"o){1 - exp[ / tr(A(/x) - ^„(/x))]}| < C\\A - A^i, 

where C is another constant depending only on the curves F and L and the operator A, 
then we have 

\DAM-DA{l)\<CKe-^\ 

Now we consider the case T[^] is not invertible. According to there is a (p with 
only finitely many non- vanishing Fourier coefficients such that T[ijj + £(j)] is invertible for 
all sufficiently small nonzero e. Let r]^ = + ecp, then the previous process shows that 

\\PnH[r],]H[rfMPnT[ve+']Pn - H[r,,]H[M]\\i < ^e""", 

where K, a are independent of e since (p has only finitely many non- vanishing Fourier 
coefficients. We consider e on the boundary of a disk D 



m&x\DA^M) - DA,e{^)\ <CKe- 
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where 



e J' 



Then the maximal modulus theorem says that 

\DAM)-DA{l)\<CKe--''. 

Since G[ip] = 1, we have 

|detr„[V'] - lim detr['0]| < Ke""", 

n— >-oo 

where < K,a < oo are constants. □ 

Lemma 5.2. Assume each fundamental domain has symmetric edge weights with respect 
to a center vertex, as illustrated in Figure 8. Assume {po,Qo) is a vertex on the boundary 
of a fundamental domain centered at (0,0). If i is even, 

Proof. Let be a layer of cylinder consisting of t fundamental domains. In other words, 
if we use the number of fundamental domains as a measure for the size of W , then the 
circumference of W is t, and the width of is 1. The transfer matrix Qt is defined 
to be a square matrix whose rows (columns) are labeled by the upper(lower) boundary 
configurations of W . Namely, let X be an upper boundary configuration, and Y be an 
lower boundary configuration, then 

Qt{X,Y) = ^ exp[-/3(^ ^ Jab(^a<^b+^ XI Jcdf^cCTd 

{(7r\reW\dW} {a,b\ar^b,(7a,crb&X} {c,d\c~d,<7c,<7deY} 

e£E{W\dW) 

Without the external magnetic field, 

{^Poqo'^Po+im,qo+jn)p = {(^poqo'^Po+imjqo+jnjp ~ {^poqo)p{'^Po+i'm,qo+jn)p 

~ ^{{PpoqoPpo+i'>'n,qo+jn)p ~ {Ppoqo)p{Ppo+im,qo+jn)p)- 

Hence it suffices to prove that 

{PpoqoPpo+im,qo+jn)p — {Ppoqo) p{Ppo+im,qo+jn) p ^ {Ppoqo Ppo+im,qo) p — {ppoqo)p{Ppo+im,go)p- 

Assume we have an s x t torus, where s and t are the number of fundamental domains 
along each direction, then 

^Ppoqo)st E,eM-PEiuMu^.}Juvaua.] HQI) ' 

where Ft is a diagonal matrix with 



Ft{X,X) = 



1 if (Jpoqo = 1 in X 
else 
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Since Qt is a symmetric matrix, all its eigenvalues are real, there exists an orthogonal matrix 

St^ such that Qt = S^^QfSt is diagonal, and the modulus of its eigenvalues decreases along 
the diagonal. Let Ft = S^^^FtSt- Moreover, since Qt has strict positive entries, the Pcrron- 
Frobenius Theorem says that the largest eigenvalue in modulus of Qt, Xi,t, is simple and 
strictly positive. Hence 

li;n(Pp.,J,,= lim^^=F.(M). 

Similarly, 

{Ppo+im,qo+jn) st = tr{QlGtQl 

and Gt is a diagonal matrix with 

Q /j^ _ I 1 ifc''po+im,qo+jn = 1 in X 

*^ ' ^ \ else 
Assume Gt = S^^GfSt, then 

lim {ppQ^ijn,qo+jn)st = Gt{^, !)• 
s— >oo 

Finally we have 

\ _ r HFtQlGtQD 

J^}^\PpoqoPpo+tm,qo+jn)st - ^^qs 

FtQjGtjl, 1) 

j:,Ft{l,k)XltGt{k,l) 

therefore 

^}^^{{PpoqoPpo+im,qo+jn)st ~ {Ppoqo) st{Ppo+im,qo+jn) st) 

Y.k^,m,m,tGt{i,k) 



Emi ^'(1' ^)Kt ^k^i Gtih k)xy ^ 



< 

= lim y {{ppQqQPpQ+im,qo)st ~ {Ppoqo)st{Ppo+im,qo)st) 

s— >-oo V 

,qo+jnPpo+im,qo+jn) st {Ppo+im,q(,+jn) st{ppo,qo+jn) st) 

= lim {ppQqQPpQ-\-im,qo)st ~ {Ppoqo) st{Ppo+im,qo) st- 
5— >oo 

The inequality follows from the Cauchy-Schwartz, it is an inner product since i is even, and 
the last equality follows from translation invariance. The result follows by letting t — t- oo, 
and the fact that the limit of the spin-spin correlation is independent of the order of s — >■ oo 
and t ^ oo. □ 
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Assume we have a finite graph. Let ^x,y be a path connecting (0, 0) and (x, y), then 
{(^GO(^xy) = Z~'^ JJ tanhJe^ (1 + — X-j^^uf^v) W (1 +tanh Je(T„cr^). 



The denominator (partition function Z) corresponds to the closed polygon configurations 
while the numerator (all the other factors on the right) corresponds to configurations on 
the graph which have an odd number of present edges at vertices (0,0) and {x,y), and 
an even number of present edges at all the other vertices. To see that, first of all, the 
expansionj^^ HeeT^ y^^'^ tanhJe """^^^ ^e^j^ ^ (1 +tanh JgCTuCrv) is the same as the weighted 
sum closed polygon configuration where all the edges except those on jxy have weight 
tanh Je, while the edges on j^y have weight ^^^^j ■ When multiplied by HeeTa: tanh Jg, 
this is the same as the weighted sum of configurations which can be obtained by a closed 
polygon configuration by changing all the present edge on jx,y to be vacant, and all the 
vacant edges on 'yx,y to be present, under the condition that all the edges in the graph 
have weight tanhj^. For all the vertices which do not lie on ■~fx,y and in the interior of 
jx,y, there are still an even number of present adjacent edges after the change. However, 
for the two endpoints (0,0) and (x, y), an odd number of adjacent edges are present in 
the new configuration. The weighted sum of new configurations can be represented by 
a dimer model on graphs by removing an vertex of the gadgets corresponding to (0, 0) 
and {x,y), where the original gadgets are illustrated as in Figure 2. Namely, the 8 local 
configurations at a degree-4 vertex, with an even number of edges present, are in one-to-one 
correspondence with dimer configurations on two gadgets, each of which is a gadget form 
the right graph of Figure 2 by removing one vertex. Figure 7 and 8 give two examples of 
such correspondence. 
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Figure 7: odd-edge configuration and modified gadget: 1 
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Figure 8: odd-edge configuration and modified gadget: 2 



Hence {(JQ.oO'x.y) can be expressed as a sum of monomer-monomer correlations the Fisher 
graph, that is, the partition function of dimer configurations on a graph with two vertices 
removed, divided by the partition function of the original graph. 

We have the following corollary 



41 



Corollary 5.3. Above the critical temperature, the monomer-monomer correlation on sym- 
metric, periodic, ferromagnetic Fisher graph decays to zero exponentially fast. The decay 
rate is rotationally invariant. 



A Toeplitz and Hankel Operators 

If is a matrix- valued function defined on the unit circle with Fourier coefficients (j)k, then 
T[0], H[(j)\ are, respectively, the semi-infinite Toeplitz and Hankel matrices, 

T[0] = 0<i,j<oo, 
HWi = 0<i,i<oo. 

(p is called the symbol of the Toeplitz and Hankel matrices. If (p is bounded, these may be 
thought as operators from the Hilbert space (.2 to ^2- In addition we writie 

~4>{z)=cl>{z-^). 

There is a simply identity relating Toeplitz and Hankel operators 

T[H] - T[4>m] = H[4>\H[i:-\ (34) 



Identity (34) is trivial. The left side has (i,j) block 

00 00 —1 00 

4>i-k-^k-j -^^(t>i-kipk-j = ^ 4>i-ki'k-j = y^^<pi+k+iip-k-j-i- 



k=0 k=-oo k=0 



Two applications of the identity ( 34 ) give (I=identity matrix) 

T[r']m = I - H[r']H[4>]. 

Let X and Y be Banach spaces. A bounded linear operator T : X ^ Y is called 
completely continuous if, for any weakly convergent sequence (x„) from X, the sequence 
Txn is norm-convergent in Y. 

A compact operator L from X to y is a linear operator such that the image under L 
of any bounded subset of X is a relatively compact subset of Y. Compact operators on 
Banach spaces are always completely continuous. 

The collection 5^p(l < p < 00) consists of all completely continuous operators A for 
which 

00 
i=i 

where Sj are singular values of A. The norm in 5p is defined by 

00 

When p = 1 the operator is called a trace-class operator. Here is an equivalent definition 
for the trace-class operator over a separable Hilbert space. A bounded linear operator A 
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over a separable Hilbert space H is said to be in the trace class if for some orthonormal 
basis {ek}^i of H, the sum 



^lli 



< {A*A)-2ek,ek > 



k 



is finite, where A* is the adjoint operator of A under the inner product in H. In this case 



is absolutely convergent and is independent of the choice of orthonormal basis. This limit 
is denoted by tr{A). When A £ ^i, this limit is equal to the sum of all eigenvalues of the 
operator A. 

A Hilbert-Schmidt operator is an operator on Hilbert space satisfying ||j4||2 < oo. A 
trace-class operator on Hilbert space is always a Hilbert-Schmidt operator. A Hilbert- 
Schmidt operator is always compact. 

Given a Hankel matrix with continuous symbol (/). The Fourier coefficients of (p decays 
exponentially fast. In this case, the Hankel operator is a Hilbert-Schmidt, and the product 
of two Hankel operators is a trace-class operator. Since r[0]r[(/)~^] = I — H[<p]H[(f)~^], we 
have T[(/)]r[(/)~^] is an operator differing from the identity by an operator of trace class, 
hence the determinant E in Widom's theorem (Lemma 4.4) is well defined. 

B Uniqueness of Gibbs Measure on Dimer Models of Non- 
bipartite Graphs 

One correspondence between the Ising spin configurations and dimer configurations is de- 
scribed in [5J. Consider the dual graph Q* of the graph Q where the original Ising model 
is defined. Replace each vertex of the dual graph by a gadget. 



Assume all the edges which do not correspond to the edges of the dual graph have weight 
1. For those edges corresponding to edges of the dual graph, assume an even number of 
such edges, surrounding each gadget, have weight less than 1, all the other dual edges have 
weights greater than or equal to 1. For example, in Figure 9, edges 61,62,63,64,65 have 
corresponding edges in the dual graph(dual edges), while the edge 60 does not correspond 
to any edge in the dual graph at all. The two-to-one correspondence between the Ising 
spin configurations and the dimer configurations is described as follows: 

1. If two adjacent spins have the same sign, and the weight of the edge of G* separating 
the two spins is bigger than 1, then the edge is present in the dimer configuration; 

2. If two adjacent spins have the same sign, and the weight of the edge of Q* separating 
the two spins is less than 1, then the edge is not present in the dimer configuration; 




k 




Figure 9: Dual Ising Graph and Decorations 
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3. If two adjacent spins have the opposite sign, and the weight of the edge of Q* separat- 
ing the two spins is less than 1, then the edge is present in the dimer configuration; 



4. If two adjacent spins have the opposite sign, and the weight of the edge of Q* sep- 
arating the two spins is bigger than 1, then the edge is not present in the dimer 
configuration; 

For example, in Figure 9, assume the weights satisfy < 1, < 1, and Wei,We3,We^ > 
1, then an Ising spin configuration on the left graph should correspond to a dimer config- 
uration on the right graph, according to the 4 criteria listed above. 

Corollary B.l. For a graph obtained from the dual graph of a ferromagnetic Ising model 
as described above, there is a unique translation-invariant Gibbs measure defined on dimers. 

Proof. It suffices to prove the result on any finite cylindrical set. Let pe be the variable 
associated to an edge e, that is, if e is present in the dimer configuration, pe = 1, otherwise 
Pe = 0. If e is not an dual edge, the configuration of e is uniquely determined by dual 
edges. For example, in Figure 

Peo = Pe3pe4(l " Pes) + Pegll " ^64)^65 

Hence it suffices to prove the result on an arbitrary cylindrical set consisting of dual edges. 
When e is a dual edge let (T„, (7„ be spins on endpoints of the dual edge e*, then 



Hence [pe-^ ■ ■ ■ Pck) depends only on spin-spin even correlation functions. The uniqueness 
of spin-spin even correlation functions implies the uniqueness of the translation-invariant 
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